GLITCHES IN THE X-RAY PULSAR IE 2259+586 
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ABSTRACT 

Starquakes are considered for fast-rotating magnetic white dwarfs. The 
X-ray pulsar IE 2259 + 586 may be such a white dwarf. It is shown that in this 
case starquakes may be responsible for the decrease of the mean spin-down rate 
Q>^ . which was observed for IE 2259 -|- 586 between 1987 and 1990. The required 

G^ I mass of the white dwarf which is identified with IE 2259 -f- 586 is ~ 1.4 — 1.5M0, 

^ . making this X-ray pulsar the most massive white dwarf known. 

o : 
Q: 

' The X-ray pulsar IE 2259+586 is located close to the center of the young 

^ ■ supernova remnant G109.1 — 1.0 whose age is estimated to be (1.2 — 1.7) x 10^ 

I I yr (Gregory and Fahlman 1980; Hughes, Harten and van den Bergh 1981). The 

^ ■ period of the pulsar, P = 6.98 s, is suitable for an ordinary binary X-ray pulsar. 

I I But, IE 2259 + 586 exhibits several peculiarities: (1) the period shows an 

^ ■ unusually stable spin-down trend with a small rate of 5.9 x 10~^^ s s~^ (Davies, 

^ ! Wood and Coe 1990); (2) no orbital Doppler modulation of the pulsar period 

• has been discovered (Koyama et al. 1989); (3) there is no firm evidence for 

'. an optical counterpart (Davies and Coe 1991), and (4) it has an unusually soft 

^ ■ X-ray spectrum. Therefore, Davies et al. (1989) have suggested that IE 2259 + 

r-| ! 586 is a single star. However, this X-ray pulsar cannot be a single neutron star 

CL|' in which, like a radio pulsar, spin-down power is a source of the radiation. The 

Q ! point is that for a neutron star with the observed spin period and the spin-down 

?-j I rate, the spin-down power cannot be larger than 10^^ ergs s~^, while the soft 

. X-ray luminosity of IE 2259 + 586 is ~ 10^^ ergs s~^ (Fahlman and Gregory 

?i ; 1981; Davies, Wood and Coe 1990). 

■ It was suggested by Paczyhski (1990) that the enigmatic X-ray pulsar IE 

'k/i \ 2259 + 586 is a single fast-rotating white dwarf with a strong magnetic field, 

' -Bg ~ 10^ Gauss, at its surface. The required mass of the white dwarf is higher 

', than 1.32M0. This lower limit is because of a centrifugal instability for the 

white dwarf if its mass is smaller than 1.32M0. Such a massive white dwarf 
may be the product of a white dwarf merger. Since the moment of inertia for a 
massive white dwarf is more than ~ 10^ times larger than for a neutron star, the 
spin-down power of the white dwarf may be as high as ~ a few xlO'^^ ergs s~^, 
which is more than enough to explain the X-ray luminosity of IE 2259 + 586. 
The soft X-ray emission of IE 2259 + 586 can be explained by the radiation 
of the hot gas in the polar caps (Usov 1993). This gas is heated by the flux of 
positrons which are created in the magnetosphere of the white dwarf. 

It was shown recently (Iwasawa, Koyama and Halpern 1992) that the spin- 
down rate, which was consistent with a constant value for the first 10 years of 
observations, from 1978 to 1987, has decreased by a factor of 2 between 1987 and 
1990. Using this, Iwasawa, Koyama and Halpern (1992) have concluded that 

the single white dwarf model of Paczyhski (1990), which required a stable P, is 
ruled out for IE 2259 + 586. In this paper I have argued that this conclusion 



is premature. The point is that if an essential part of the white dwarf is sohd, 
starquakes have to have happened occasionaUy in the process of deceleration 
of the white dwarf rotation. In turn, these starquakes may result in sudden 
changes of P, which are similar to the pulsar glitches. One of these changes of 
P may be responsible for the decrease of the mean spin-down rate which was 
observed between 1987 and 1990. 



2. STARQUAKES AND THE PULSE PERIOD EVOLUTION 

Let us consider a fast-rotating magnetic white dwarf with a solid core. 
As the rotation of the white dwarf slows down because of the electromagnetic 
torque, centrifugal forces on the core decrease, and gravity pulls it toward a less 
oblate shape, thereby stressing it. When stresses in the core reach a critical 
value, the core cracks, some stress is relieved, and the excess oblateness, due 
to the core rigidity, is reduced. The moment of inertia of the white dwarf 
is suddenly decreased, and by conservation of angular momentum the angular 
velocity of the white dwarf rotation, Q = 27r/P, is suddenly increased. This 
is similar to the starquakes which were considered for neutron stars to explain 
the pulsar glitches (Ruderman 1969; Baym and Pines 1971; Pines, Shaham and 
Ruderman 1972). 

The oblateness of a rotating white dwarf may be characterized by the 
following dimensionless parameter: 

(1) 







where / is the moment of inertia of the white dwarf, and 1^ is its moment of 
inertia in the case that the white dwarf does not rotate. 

Since the angular momentum of the white dwarf is conserved in the quake, 
the sudden change in the oblateness is 

A£=^ = ^' (2) 

where A/ is the change of /, and AP is the jump of P. 

The characteristic time between quakes is (Baym and Pines 1971) 

h^^s^\^^\^ (3) 

where = P/P is the time that characterized the rate at which the white dwarf 
slows down due to loss of rotational energy [r^ is ~ 3 x 10^ yr for IE 2259 + 

586 ], 

= pc/Anip is the ion density, pc is the density of matter in the white dwarf 
core, A is the atomic weight of ion, Z is its electric charge, rup is the mass of 



proton, e is the charge of electron, G is the gravitation constant, Mc is the mass 
of the white dwarf core, and i?c is its radius. Here and below we assume that 
the main part of the white dwarf mass is solid. To get equations (4) and (5) 
for D and B the solid core was considered as a rotating self-gravitating sphere 
composed of incompressible matter of uniform shear modulus. The estimate of 
tg by means of equation (3) is uncertain within a factor of 2-3 or so. 

The more reasonable parameters of the white dwarf which is identified with 
IE 2259 + 586 are the mass of the white dwarf M ~ 1.45M0 (see below), the 
stellar radius ~ 3x10^ cm, ~ IM©, ~ 2x10^ cm, = Mc/(47r/3)i?3 ~ 
6 X 10^ g cm~^, and 1^ ~ 10^° g cm^. Substituting these parameters into 
equations (2) - (5), for the white dwarf which consists of iron, Z = 2Q and 
A = 56, we obtain D ~ 10^*^ ergs, ~ 7 x 10^^ cm"*^, a ~ 1.4 x 10~^° cm, 
B ~ 10^^ ergs, u;^ ~ 20 s~^, tg ~ 7 x 10^ | Ae \ yr. Taking into account equation 
(2), we have the following relation between tq and | AP \ /P for IE 2259 -|- 586: 

I AP I 

tg~7xl0^^^yr (6) 

The available observational data on the changes of the pulse period of IE 
2259 -|- 586 from 1978 to 1990 are given in Figure 1. This data may be in 
agreement with the single white dwarf model of IE 2259 -|- 586 if we assume the 
existence of glitches with AP = — (1 — 2) x 10~^ s and ~ 3 — 10 yr or more 
(see Fig. 1). For this value of AP equation (6) yields tg ~ 10 — 20 yr, which 
might be uncertain within a factor of 2-3. This value of tq is consistent with the 
data on the evolution of the pulse period of IE 2259 -|- 586. 

One of the main assumptions which we have used in this paper is the 
existence of massive solid core inside the white dwarf which is identified with IE 
2259 -|- 586. The cooling of white dwarfs and their crystallization were considered 
mainly for white dwarfs with mass M ~ O.6M0, the "typical" mass remnant of 
single-star evolution (for a review, see D'Antona and Mazzitelli 1990). For such 
a white dwarf the time, Tcr, from its formation to the stage of crystallization is 
~ 10^ — 10^ yr. This is much more than the age of the X-ray pulsar IE 2259 -|- 
586 in the model of Paczy'nski (1990) in which the pulsar cannot be much older 

than Ts = P/P ~ 3 x 10^ yr, but it may be as young as ~ 1.5 x 10^ yr, if it 
is related to the supernova remnant G109.1 - 1.0. The value of drops with 
increasing M, and it may be as small as ~ 10^ yr for the white dwarf with the 
mass M ~ 1.4 — I.5M0 which relates to IE 2259 -|- 586. Indeed, for a plasma 
containing only one species of ion, the temperature of crystallization is (see, for 
example. Van Horn 1971; D'Antona and Mazzitelli 1990) 

T„ = 2.28xl0'r-|l(^j^^)"' K, (7) 

where F is the dimensionless parameter which is somewhere between 64 and 210. 
The more reasonable value of F is 150. 

Substituting F = 150, Z = 26 and ^ = 56 into equation (7) we have 



The temperature inside a massive white dwarf with the age ~ 10^ yr is not 
higher than ~ 10^ K (Vila 1966; Savedoff, Van Horn and Vila 1969). A pure 
iron matter with this temperature has to crystallize if its density is equal or 
more than ~ 5 x 10^ g cm~^ (see equation (8)). Both the mean density and 
the density at the center of a white dwarf increase very sharply when the white 
dwarf mass goes to the edge of stability (Geroyannis and Hadjopoulos 1989 and 
references therein). Using the paper of Geroyannis and Hadjopoulos (1989), we 
can see that the main part of the white dwarf mass is inside of the surface at 
which the density is equal to 5 x 10^ g cm~^ only if the mass of the white dwarf is 
~ 1.4- 1.5M0. Therefore, if the white dwarf, which is identified with IE 2259 + 
586, consists of iron and its mass is near the edge of stability, M ^ 1.4 — 1.5M0, 
the assumption that the main part of the white dwarf mass is solid is reasonable. 

3. DISCUSSION 

In this paper I have suggested the existence of glitches with the amplitude 
AP ~ —(1 — 2) X 10~^ s and the characteristic time between glitches ~ a few 
x(l — 10) yr for the enigmatic X-ray pulsar IE 2259 + 586. This is in agreement 
with the single white dwarf model of Paczyhski (1990) for IE 2259 + 586. The 
observation of glitches would be a confirmation of this model. 

The other way to verify the single white dwarf model is the observation of 
7-rays from IE 2259 + 586. The point is that fast-rotating white dwarfs, P ~ 
a few s, with strong magnetic field, ~ 10^ — 10^ Gauss, are similar to the 
radio pulsars, and both the particle acceleration to ultrarelativistic energies and 
generation of 7-rays have to be in their magnetospheres (Usov 1988, 1993). If 
IE 2259 -|- 586 is such a white dwarf, the expected fiux of 7-rays with energies 
~ 10^ - 10^ MeV may be high enough to be detected by EGRET (Usov 1993). 

FIGURE CAPTION 

Figure 1. Pulse period histiry of IE 2259 586 from 1978 to 1990. The 
solid line is the best fit to the points before 1989 without glitches. The dashed 
lines show the fit to the points with glitches. The amplitude of the glitches is 
somewhere between 10~^ s and 2 x 10~^ s. The characteristic time between 
glitches may be as small as a few years (see Fig. lb). In this case the spin-down 
rate between glitches is essentially more than the mean spin-down rate for the 
first 10 years of observations. 
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